Entanglement of an impurity and conduction spins in the Kondo model 
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Based on Yosida's ground state of the single-impurity Kondo Hamiltonian, we study three kinds 
of entanglement between an impurity and conduction electron spins. First, it is shown that the 
impurity spin is maximally entangled with all the conduction electrons. Second, a two-spin density 
matrix of the impurity spin and one conduction electron spin is given by a Werner state. We find 
that the impurity spin is not entangled with one conduction electron spin even within the Kondo 
screening length although there is the spin-spin correlation between them. Third, we show the 
density matrix of two conduction electron spins is nearly same to that of a free electron gas. The 
single impurity does not change the entanglement structure of the conduction electrons in contrast 
to the dramatic change in electrical resistance. 
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Entanglement, the quantum correlation between sub- 
systems, is considered to be one of the key concepts 
in quantum mechanics and quantum information sci- 
ence 0. A study of entangled structures of quantum 
many-body systems is of great importance for provid- 
ing us not only a new view on their physical properties 
but also the basic knowledge for fabricating quantum in- 
formation processors. For example, entanglement of the 
ground states of one-dimensional quantum spin lattice 
models has been intensively investigated in context of 
quantum phase transition 0,0. In a non- interacting 
electron gas, entanglement of two electron spins due to 
the Pauli exclusion principle has been studied in connec- 
tion with its characteristic length scale, the Fermi wave- 
length Xp 0,01 El- Entanglement of two electron spins of 
a BCS superconductor has also examined with relation 
to the coherence length £ 0- 



The Kondo model describing the exchange interac- 
tion between the impurity spin and the conduction elec- 
trons has been one of challenging quantum many-body 
problems in condensed matter physics 

00 EH Vari- 
ous theoretical tools such as Anderson's scaling theory, 
Wilson's numerical renormalization group approach, and 
the Bethe ansatz, etc. have been applied to the Kondo 
model 0. Recently, Kondo effects in nanodevices have 
been revisited with potential applications to spin de- 
vices or spin qubits 0, 0, [nj. Although the Kondo 
model have been very intensively studied, its entangle- 
ment structure remains to be explored. In this paper 
we investigate three kinds of entanglement between the 
impurity spin and conduction electron spins as shown 
in Fig. ^ From this study on the entanglement struc- 
ture of the Kondo model we provide a clear view on 
the Kondo screening cloud with the size of the Kondo 
screening length which is the holy grail in the Kondo 
physics 0, ■ Also we discuss a possibility for the use 
of an electron gas in coupling spin qubits [lfj, [ij . 
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FIG. 1: (color online). Schematic of three kinds of entan- 
glement: entanglement (a) between the impurity spin and all 
the conduction electrons, (b) between the impurity spin and 
one conduction electron spin at r, and (c) between two con- 
duction electrons at n and V2- The rectangular box is the 
conduction electrons, the arrow with a filled circle represents 
the impurity spin, and the red and blue arrows denote the 
conduction electron spins. 



Let us consider the spin-i Kondo Hamiltonian 



H = J2 - 2 JS • s(0) . 
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where S is the spin of the magnetic impurity, s(0) the spin 
of the conduction electrons at r = 0, and cj^ creates a 
conduction electron with momentum k and spin a. The 
exchange interaction J is assumed to be negative, so the 
impurity spin is anti-parallel to the conduction spins. 

Yosida presented the variational gro und state of Eq. Q 
given by a singlet state 0, Q, [H| 
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where |x|) is the spin-up state of the impurity spin and 
\4>l) denotes the state of conduction electrons with one 
excess down spin. The simplest form of \<f>^) is the state 
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with one down spin added outside the filled Fermi sphere 



\4>i) 
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fe>fej; 



where are variational parameters, J\f the normaliza- 
tion factor, and \F) — Y\ k<k c^JO) the filled Fermi 
sphere. The variational parameters Tk are given by 



(4) 



where E B = k B T K = D exp[2/3JJV(0)]. Here T K is 
called the Kondo temperature, 2D the bandwidth of the 
conduction electrons, and N(0) the density of states of 
the conduction electrons. The characteristic length scale 
of the Kondo model is the Kondo screening length £r-, 
the size of the Kondo screening cloud 
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where vf is the Fermi velocity, Ep the Fermi energy, 
and k F the Fermi momentum. Since Ef ^> k B Tx, S,k is 
much larger than the Fermi wavelength Xf = 2n/kF- 

(i) Entanglement between the impurity spin and the 
conduction electrons. As depicted in Fig. [fl(a), let us 
investigate entanglement between the impurity spin and 
all the conduction electrons. Since the total state, Eq. @ 
is a pure state, entanglement between them is quantified 
by the von Neumann entropy of the reduced density ma- 
trix p im of the impurity spin, S(p im ) = -Tr [pi m logp im ]. 
Here the subscript 'im' represents the impurity. By trac- 
ing out the degrees of freedom of the conduction electrons 
in Eq. (J2J), one easily obtains p lm = Tr con | v E's)( v I's|. It is 
given by the fully mixed state 
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which gives the maximum entropy S{p- lm ) = 1. Thus we 
find that the impurity spin is maximally entangled with 
the conduction electrons. 

(ii) Entanglement of the impurity spin and one con- 
duction electron. If a magnetic impurity is immersed in 
a Fermi sea, conduction electrons will align in order to 
screen the magnetic field by the impurity. To put it sim- 
ply, this magnetic disturbance is called the Kondo screen- 
ing cloud of electrons (or an electron) with the size of 
order forming a singlet with an impurity (l2^. So one 
may expect some correlation or entanglement between 
the impurity spin and a conduction electron in the Kondo 
screening cloud. 

Let us consider the impurity spin and one conduction 
spin at r from the impurity as shown in Fig. E(b). En- 
tanglement between them can be measured if one obtains 
the two-spin density matrix by tracing out all the degrees 
of freedom of Eq. (J2J except the impurity spin and one 



conduction spin at r. In the second quantization, it is 
given as follows 

pSo'/J'W - 5(*.l^(')^(O)0a(O)^(r)|*,) , (7) 

where an operator ^t(r) creates a conduction electron 

with spin (3 at r and <pa is a creation operator of the 
impurity spin a at origin. Here a, a',/3, [3 1 refer to the 
spin indices, i.e., j or |. 

From Eqs. (J2J and Q, we obtain 



P (2) (r)= 1 - 



|0fT;tT 

PiV-M P1T;1T 




(8) 







where the diagonal elements are given by 

PTT;TT = (^iI^tW^tWI^) - 
PU;U = ( < ^llV'I( r )V'l( r )l0i> » 
PIV.I1 = (<HV^( r )V'T( r )l ( / , T) . 

Piun = (<H^I( r )^i( r )l<?H) . 

and the off-diagonal elements are given by 

PU;U = -<0il$(r)lM r )l0T> = P*VAl 
Since the state |</>|) has N/2 up and N/2 + 1 down spins 
of the conduction electrons, the matrix element /OfTiTT ^ s 
given by the electron density with up spin 



(9a) 
(9b) 
(9c) 
(9d) 

(9e) 



N n 
PTT;TT = W " j 



(10) 



where V is the volume of the conduction electrons and n 
is the electron density. Similarly we have = n /2 ■ 

The matrix element PtX;U ^ s given by 



Pn-,n 

where f(r) is defined by 
1 



VN . 
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(12) 



kk'>fcF 



The system is isotropic, so /(r) depends only on the dis- 
tance r. We get the off-diagonal element PfJ-iiT = —f(f) ■ 
Thus is given by the sum of the fully mixed state 
and a spin-singlet state 



P {2) = 







f+/(r) 
-f(r) 








f + /(r) 
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= n-+/(r)|^-))(^-)| 
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where |*( - )) = 4=(\1 |) - || t)) and I is the 4x4 unit 
matrix. The meaning of Eq. (|13|l is that the first term 
is the uniform back ground and the second is the singlet 
state forming a Kondo screening cloud. 

Since Trp( 2 ) = n+/(r), let us introduce the normalized 
two-spin density matrix p defined by p^ = (n + f)p ■ 
Then we find that p is given by the Werner state 0, 



(l-p)-+p|*(-))(*(-)| 



(14) 



where p = f j(n + /). The impurity spin is entangled 
with one conduction electron spin at r if p > 1/3, that 
is, / > n/2. 

Let us examine whether p is entangled or not by calcu- 
lating /(r). Since r k is real, we rewrite /(r) = Vf(r) 2 /JV 
where the function /(r) is defined by 



/» 4^ rkei 



(15) 



After some calculations, f(r) becomes the integral form 
• D/E F sin 



N(0) 
kpv 



k F ry/l + tE B /E F 



t + 1 



dt, (16) 



where A^(0) is the density of state at Fermi level and be- 
comes iV(0) = rnkp /2ir 2 h 2 ~ 3n/4E F for a free electron 
gas. The normalization factor Af is given by 



k>kp 



VN(0) 
E B 



y(D,E B ,E F ) 



(17a) 



where the definite integral y(D, Eb, Ep) is defined by 

"D/ Ef 



y(D,E B ,E F ) = 



Thus we obtain 



f(r) = 



3n E B 



^l + tE B /E F 
f N {r) 2 



4 E F y{D,E B ,E F ) 



dt . (17b) 



(18) 



where /jv(r) = f(r)/N(0). The condition of entangle- 
ment of p, f > n/2, becomes 



fN{rf 



2\E F ) y{D,E B ,E F ) 



> 1 . 



(19) 



For usual Kondo systems, the Kondo temperature Tk 
ranges from few to hundreds Kelvin, Tk — 1 ~ 300 K. 
The Fermi temperature Tp = Ep/ks is about Tp ~ 
10 4 K. This implies E B /E F ~ 10~ 2 - 10~ 4 . How- 
ever, according to Eq. JJJ, the small ratio of Eb/E f = 
^bTkIEf indicates the Kondo screening length £k is 
much larger than the Fermi wavelength . If the band- 
width D is assumed to be Eb <C D -^i E F , we have 
y(D, E B , Ep) ~ 1. Also /jv(r) ~ 0{l) as shown in 
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FIG. 2: /jv(r)//jv(0) as a function of kpr. The numerical 
values, Ef, Eb and D are taken in the text. 



Fig. (0. Then f/n could not be greater than the half, 
that is, Eq. i|19|) is not satisfied. Therefore the impurity 
spin is not entangled with one conduction electron spin 
even within the Kondo screening length, i.e., r < ^k- 

Fig. Q shows the numerical calculation of 
/jv(r)//Ar(0) as a function of k F r. /jvW extends 
over the Kondo screening length £k- We find that f(r) 
is noting but the spin-spin correlation function between 
the impurity spin and the conduction spin at r, 

(fini^M) = PTT;TT + PlUll ~ PU;T1 _ PlT;iT 

= -2/(r). (20) 

Notice that f(r) is qualitatively similar to the spatial 
spin-spin correlation function (S • sfr)) = | (crf m (T z (r)), 
which was calculated by Chen et a l 1 161 and by Guber- 
natis et al for the Anderson model [17j . It is interesting 
that p( 2 Hs given by a pseudo entangled state in liquid- 
state NMR [H G]J. The Kondo screening cloud [aH^ 
could be detected through the extra Knight shift experi- 
ment which measures f(r). 

Our result shows that even if there is the spin-spin 
correlation between the impurity and a conduction elec- 
tron in the Kondo screening cloud, entanglement between 
them vanishes. A simple explanation for this is as follows. 
Since many conduction electrons are coupled to a single 
impurity, the quantum correlation between the impurity 
and each conduction electron is very tiny so there is no 
entanglement between them. 

(Hi) Entanglement between two conduction electron 
spins. It has been believed that conduction electrons 
within the Kondo screening cloud are mutually corre- 
lated because they have information on the same impu- 
rity |9j. We examine whether the impurity spin induces 
the non-classical correlation, i.e., entanglement between 
conduction electrons. 

Consider two conduction spins at r*i and r2 as shown in 
FigH-(c). The density matrix p con of all the conduction 
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electrons is given by the mixture of and \(j>±) 

Peon = Tr im |* 5 )(* s | = 1 ( 10 T ><0 T | + ) . 



(21) 



From Eq. (|21[1 , the density matrix of two conduction 
electron spins is written by 



~(2) / \ 

" a x ,<r 2 \a\ ,cJ2 v ' y 



i T r 

2 



^( r 2)V'i' 1 ( r l)^l( r l)^2( r 2)Pco 1 

After a lengthy calculation, we obtain 

p( 2 )(r 1 ,r 2 )=^ e (r') + Ap(r 1 ,r 2 ) t 



(22) 



(23) 



where r = ri — r 2 and Pf r j e (r) is the density matrix of 
two electron spins of a free electron gas 0, 



1-5 2 

1 -g 2 

-g 2 1 

1-g 2 



(24) 



with g(r') = -fr J2k<k e • The smau change is 



A/?(ri,r 2 ) 



n 2 3E B 
~8~ 2Ef~ 



a + b 

a b 

6 a 

a + 6 



(25) 



where a = ,g(0) [ f N {r x ) 2 + f N (r 2 ) 2 ] /2 and & = 
.9( r ') /W( r i) ijv( r 2)- Since Eb/ Ep is very small, we have 
p (2) (ri,r 2 ) « p| r 2 ) (ri,r 2 ). Of course, if r 1 ,r 2 > £k, then 
one gets p (2) (ri, r 2 ) = p[ r 2 c } c (ri, r 2 ). 

This result implies that the entanglement structure of 
the conduction electrons is little affected by the single 
impurity. In contrast to the traditional belief Q , the mu- 
tual correlation between conduction electrons is mainly 
due to the Pauli exclusion principle not due to the impu- 
rity. So their correlation length is the Fermi wavelength. 
The reason is that the number of conduction electrons 
are very large so they act as a reservoir. Although |<^>j) 
is orthogonal to \(f>i), \4>t) differs from |<^>j_) and from \F) 
by a single spin. 

Before conclusion, we would like to mention an inter- 
esting experiment which reported the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction of two spin qubits on 
quantum dots mediated by an electron gas [10l 1 1 1 1 | . The 
RKKY interaction may be used to couple spin qubits. 
However, our result implies that a spin qubit could be 
maximally entangled with the electron gas and become 
fully mixed if one does not deal with a quantum sate of 
a spin qubit and the electron gas as a whole. Our an- 
other study |l9| on a two-spin boson model, similar to 
the two-impurity Kondo model, suggests that the elec- 
tron gas could be used for an interaction mediator if the 
electron gas are separable from two spin qubits. 



Our analysis is based on Yosida's ground state which 
describes the Kondo physics in an simple but effective 
way. There are another ways to solve the problem: 
the variational ground state of the Anderson Hamilto- 
nian [20I El| , the quantum Monte Carlo method ^3 1 the 
perturbative renormalization |l6l| , and the density matrix 
renormalization group for a tight-binding model |l2j |. 

In conclusion, we have studied entanglement of the 
impurity spin and the conduction electron spins in the 
Kondo model based on Yosida's ground state. First, it 
has been shown that the impurity spin is maximally en- 
tangled with all the conduction electrons. Second, the 
two-spin density matrix of the impurity spin and one con- 
duction electron spin at r is given by a Werner state. It 
has been found that the impurity spin is not entangled 
with a conduction electron spin within the Kondo screen- 
ing cloud even though there exists the spin-spin correla- 
tion between them. Third, the entanglement structure of 
the conduction electrons is little affected by the impurity 
in contrast to the strong effect on electrical resistance. 
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